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1 Introduction 

The study of g-factors of subatomic particles can trace its roots back to the 1921 
paper by Stern B which was a proposal to study space quantization with an apparatus 
which is now called a "Stern-Gerlach apparatus". By 1924 the famous experiments 
had been done|| and a review paper was written summarizing their results. Their 
final conclusion, that "to within 10% the magnetic moment of the electron was one 
Bohr magneton", meant in modern language that the g- value of the electron was 2, 
where the gyromagnetic ratio g is the proportionality constant between the magnetic 
moment and the spin, 

The discovery that g e ^ 2,|| and the calculation by Schwinger|6[] predicting that 
(to first order) the radiative correction to g e was cx/tt, were important early steps in 
the development of Quantum Electrodynamics (see Fig. [I]). 





Figure 1: The Feynman graphs for g = 2 and the lowest order radiative correction 
(Schwinger term) 



The long lifetime of the muon permits a precision measurement of its anomalous 
moment at the ppm level. The muon magnetic moment is given by|7] 

eh (g — 2) 

fXf, = (1 + a n)-^- where a M = — - — . (2) 
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The electron anomalous magnetic moment has been measured to a few parts 
per billion, || which can be completely described by QED of electrons and photons 
to eighth order, (- ) 4 . The contributions of virtual muons, tauons, etc. enter at 
the few ppb level. The calculation of the electron anomalous moment is limited by 
the knowledge of the fine-structure constant. With the reliability of modern QED 
calculations, Kinoshita has turned things around and has used the electron g value 
measurement to give the best value for a.|| 

The relative contribution of heavier particles to the muon anomaly scales as 
(m M /m e ) 2 and the famous CERN experiment, |i0| which obtained a relative error 



on a M of ±7.3 parts per million (ppm), easily observed the predicted ~ 60 ppm con- 
tribution of virtual hadrons. 

In 1984 efforts began to make a new measurement of the muon anomalous moment 
to a precision of ±0.35 ppm, which would represent a 5 standard deviation observation 
of the electroweak contribution, and would also be sensitive to contributions from 
"new physics" such as muon substructure or supersymmetry. 



2 Theoretical Contributions to (g — 2) 



The standard model value of a M is given by a M (SM) = a^(QED) ± ^(hadronic) ± 
a M (weak) and any contribution from new physics would be reflected in a measured 
value which did not agree with the standard model. 

Comparison of the measurements and calculations of the electron g value gives 
one great confidence in our understanding of QED0 to the level needed for muon 
(g — 2). Taking the value of a from the electron (g — 2),|| yields the total QED 
contribution a M (QED) = 116 584 705.7(1.8) (0.5) x 10- n [|II 

The hadronic contribution to (g — 2) cannot be calculated directly, but must 
be determined from data. The first-order hadronic vacuum polarization dominates 
the uncertainty in the theoretical value of a M , since it is calculated using disper- 
sion theory and data from e + e _ — > hadrons and hadronic r decay as input. The 
various order hadronic contributions are shown in Fig. |2]. Diagrams for hadropro- 
duction and hadronic r decay are shown in Fig. ^. The most precise determination 
of the first-order hadronic contribution is [12] a^(had; 1) = 6924(62) x 10~ n which is 
59.39 ± 0.53 ppm of a M , but there is continuing discussion of the use of CVC and the 

-101(6) x 10~ n 



.(had; 2) 



r-decay data. [13) The higher-order contribution is|14j 

The hadronic light-by-light scattering shown in Fig. f| has now been calculated 
by two groups, [15], [Uj using essentially the same model, and agreement is found: 

However the two groups disagree on the uncertainty on 

The uncertainty 



t (had;lbl) = -85(32)10 



-ii 



16 



the calculation, and I have taken the larger error from Ref. 
in this contribution could be reduced substantially by the appropriate calculation on 
the lattice, and perhaps by other additional calculations as well.[|T7| 
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Figure 2: (a) The lowest-order hadronic contribution, (b-d) Higher order hadronic 
contributions except for the light-by-light scattering contribution. 





Figure 3: (a). The hadroproduction process which enters the dispersion relation, (b) 
Hadronic r decay. 



The total hadronic contribution is given by a^(had; 1 + 2 + lbl) = 6738(70) x 10 -11 
which is 57.79 ± 0.60 ppm of a M , with an uncertainty dominated by the uncertainty 
on the first-order hadronic vacuum polarization. It is precisely this contribution 
which is being addressed by the programs to measure R(s) at BES and the Budker 
Institute. [|T^] We look forward to additional high quality data from these experiments, 
from DAPHNE, 
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as well as r-decay data from CLEO to further reduce the uncer- 
tainty on the hadronic contribution. 

The standard model electroweak contribution arises from the diagrams shown in 




-85 (32) x 10 



Figure 4: (a). The hadronic light-by-light scattering contribution. 
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Fig. | (the standard model Higgs contribution is negligible). The single-loop W and 
Z contributions were calculated by a number of authors shortly after the standard 

11 or 1.7 ppm 



model was developed. |20], |2l 
of a 



21 The result is a M (weak; 1) = 195 x 10" 



Partial calculations p4|, |25 | of the two-loop electroweak contributions indicated 
that they might not be small. The full calculation p6|, |2Tf which was later confirmed 
independently |28j showed that the total first and second-order weak contribution was 



20% less than the first order result, 
which is 1.30 ± 0.03 ppm of a M . 



The result is a M (weak; 1 + 2) = 151(4) x 10 



-n 




Figure 5: (a,b). The single-loop electroweak contributions. (The standard Higgs con- 
tribution is negligible.) (c-e) Examples of the higher order electroweak contributions. 
(SeeRefs. EH, E7B.) 



The standard model prediction for is a M (SM) = (116 591 594.7 ± 70) x 10~ n 
(±0.60 ppm). 

A great deal has been written about the possible contribution to the muon (g — 2) 
value from non-standard model physics. Just as the proton substructure produces 
a g- value which is not equal to two, muon substructure would also contribute to 
the anomalous moment, the critical issue being the scale of the substructure. p9[ A 
standard model value for (g — 2) at the 0.35 ppm level would restrict the substructure 
scale to around 5 TeV. 

In Fig. |5](a) the triple gauge vertex WW'j appears, and it is through this diagram 
that the muon (g — 2) obtains its sensitivity to W substructure and anomalous gauge 
couplings. The combined sensitivity of LEP1, LEP2 and (g — 2), and the unique 
contribution which (g — 2) makes in constraining the existence of such couplings, is 
described by Renard et al. [ftt| 

Supersymmetry has become a serious candidate for physics beyond the standard 
model. The SUSY contribution is shown in Fig. H In the case of large tan j3, the 
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Figure 6: The lowest order supersymmetric contributions to (g — 2). 



chargino diagram dominates and the contribution to (g— 2) from SUSY is given by|3T 

aJSUSY) ~ %— tan/3 ~ 140 x KT 11 ( 100 GeV )2 t an/3. ( 3 ) 

8ir sin 6w m 2 m 

The goal of E821 is to reach a precision of ±40 x 10~ n (±0.35 ppm), so the factor of 
140 above corresponds to 1.2 ppm. For m = 750 GeV and tan/? = 40, a M (SUSY) = 

ioo x nr 11 . 



3 The New (g — 2) Experiment 



For polarized muons moving in a uniform magnetic field B which is perpendicular to 
the muon spin direction and to the plane of the orbit, and with an electric quadrupole 
field E for vertical focusing, 0] the difference angular frequency, u a , between the spin 
precession frequency uj s and the cyclotron frequency u c , is given by 



w fl = 

m 



(3 x E 



(4) 



The dependence of u a on the electric field is eliminated by storing muons with the 
"magic" 7^=29.3, which corresponds to a muon momentum p M = 3.09 GeV/c. Hence 
measurement of uo a and of B determines a^. At the magic gamma, the muon lifetime 
is 7r = 64.4 /xs, the (g — 2) precession period is 4.37 /is, and for the central orbit 
radius of 7.11 m the cyclotron period is 149 ns. 

The storage ring magnet is a superferric 700 ton, 14 m diameter circular "C"- 
magnet, with the opening facing inward towards the ring center. The field is excited 
by three 14 m diameter superconducting coils which carry 5.2 kA from a low voltage 
power supply to produce the 1.45 T magnetic field. [32| The short term field stability 
over several AGS cycles is better than 0.1 ppm. 

The magnetic field which enters in Eq. [| is the average field seen by the muon 
distribution. Since direct injection of muons does not uniformly fill the phase space, 



we used a tracking code to calculate the distribution of muons in the storage ring. 
The radial distribution obtained from this tracking code was compared with the 
distribution obtained from observing the beam debunching in the ring at early times. 
The two distributions agreed quite well. In the 1999 run, a straw-tube array was 
operational at one detector location which provided information on the decay positron 
trajectories coming out of the storage region. These data will permit us to reconstruct 
directly the muon spatial distribution in one section of the ring. 

In 1998, direct muon injection into the storage ring was employed for the first time. 
The AGS performance, the beamline and the inflector magnet were as described in 
Carey et al.,|33| and except for the muon injection many of the experimental details 



are the same as described there. The positive muon beam with the magic momentum 
is formed by collecting the highest energy muons from pion decay in a 72 m long 
decay section of our beamline, which results in a muon polarization of 96%. The flux 
incident on the inflector magnet was 2 x 10 6 per fill of the ring. 

The 10 mrad kick needed to put the muon beam onto a stable orbit was achieved 
with a currents, since usual magnetic kicker techniques would spoil the precision mag- 
netic field. Three pulse- forming networks powered three identical 1.7 m long kicker 
sections consisting of parallel plates on either side of the beam. Current flowed down 
one side crossed over and flowed back up the other side. The kicker plate geometry 
and composition was chosen to minimize eddy currents, and the eddy current effect 
on the total field seen by the muons was less than 0.1 ppm 20 /xs after injection. 
The current pulse, which was formed by an under-damped LCR circuit, had a peak 
current of 4100 A and a pulse base width of 400 ns. Since the cyclotron period of 
the muon beam from the AGS was 149 ns, the beam was kicked several times before 
the kicker pulse died out. With muon injection, the number of detected positrons 
per hour was increased by an order of magnitude over pion injection. Thus the use 
of a muon beam an order of magnitude less intense than the pion beam resulted in 
a substantial increase in stored muons per fill of the ring, with the injection related 
background reduced by about a factor of 50. 

Positrons from the in-flight decay fi + — > e + i / e V fl are detected with Pb-scintillating 
fiber calorimeters placed symmetrically at 24 positions around the inside of the storage 
ring|J3). The decay positron time spectrum isflOl, |3f 



N e~ thT [1 + A(E) cos (u a t + <p{E))} . (5) 

The normalization constant N and the parity violating asymmetry parameter A{E) 
depend on the energy threshold placed on the positrons. The fractional statistical 
error on oj a is proportional to A^N^ 1 ' 2 , where N e is the number of decay positrons 
detected above some energy threshold. For an energy threshold of 1.8 GeV, we 
measure A to be 0.34 consistent with its theoretical value, which we attribute to 
the good calorimeter energy resolution (a/E = 10% at 1 GeV) and a scalloped vacuum 
chamber which minimizes pre-showering before the positrons reach the calorimeters. 
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The photomultiplier tubes of the calorimeter were gated off before injection and 
when gated on, they recovered to 90% pulse height in < 400 ns, and reached full 
operating gain in several [is. With the reduced flash following injection it was possible 
to begin counting as soon as 5 fis after injection in the region of the ring 180° around 
from the injection point. 

The calorimeter pulses were continuously sampled by custom 400 MHz wave- 
form digitizers (WFDs), which provided both timing and energy information for the 
positrons. Both the NMR and WFD clocks were phase-locked to the same LORAN-C 
frequency signal. The waveforms were zero-suppressed, and stored in memory in the 
WFD until the end of the AGS cycle. Between AGS acceleration cycles the WFD 
data were written to tape for off-line analysis, as were the calorimeter calibration 
data and the magnetic field data. 

A laser /LED calibration system was used to monitor calorimeter time and gain 
shifts during the data-collection period. Early-to-late timing shifts over the first 200 
ius were on-average less than 20 ps, which is needed to keep systematic timing errors 
smaller than 0.1 ppm. 

For the offline analysis, the detector response (waveform shape) to positrons was 
determined from our data for each calorimeter. These shapes were then fit to all 
pulses in the data to determine a time, an amplitude and a width parameter for each 
pulse. Time histograms were formed for each detector. These independent data sets 
were analyzed separately and were in ag reement (x 2 /v = 17.2/20). 

A completely blind analysis was performed on the data. Arbitrary offsets were 
put on the muon frequency and the proton frequency from the NMR probes. Each 
offset was known by one person, making it impossible to determine the actual value 
of a M . Only when the analysis of both the magnetic field and u a were completed were 
the offsets removed and the new value of a M determined. 

After the offsets were removed, it was necessary to make two corrections to the 
frequency obtained from the fitting. For muons with the "magic" momentum, u a is 
not affected by the electric field. For the ensemble of muons in our storage ring there is 
a small electric field correction to u> a since not all muons are at the magic momentum. 



There is also a pitch correction because of the vertical betatron oscillations [JR], p6 
The sum of these two corrections for these data is (0.9 ± 0.2) ppm. 



The dominant systematic errors which were reported in our first measurement ||33|| 
have been completely eliminated. The remaining systematic errors are under study. 
With many of them approximated by upper limits, one obtains a total systematic 
error of < 1 ppm, with the systematic error assigned to the magnetic field of 0.5 
ppm. Since the study of systematic errors is a source of much continuing work, we 
have chosen not to present a detailed list at this time. 

We do wish to note the substantial improvement in the magnetic field quality 
which has been obtained by additional shimming. In Fig. [7| we show the average 
magnetic field from the 1997 and 1998 runs. The field uniformity over the storage 
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aperture in 1998 was almost an order of magnitude better than was obtained in the 
CERN experiment. |flO[ 
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Figure 7: The magnetic field profile averaged over azimuth for July 1997 when the 



data reported in Carey et al.[33] were taken, and for August 1998 when the new data 
reported here were taken. Each contour represents a 2 ppm change. 



4 New Results 



One month after the Symposium, we finished our analysis of the 1998 data and 
obtained a new result at the precision of ±5 ppm. 

Our experiment measures the frequency ratio R = uj a /u p , where uj p is the free pro- 
ton NMR frequency in our magnetic field. Including the pitch and electric field cor- 
rections, we obtain R = 3.707 201(19) x 10~ 3 , where the 5 ppm error includes a 1 ppm 
systematic error estimate. We obtain from = R/ (X — R) = 116 591 91(59) x 



10 10 in which A = \i v 



3.183 345 39(10). 0, This new result is in good agree- 

| and our previous 



ment with the mean of the CERN measurements for a M + and a M 
measurement of a. 



p3fl which are tabulated below. 
Assuming CPT symmetry, the weighted mean of the four measurements gives a 
new world average of a M = 116 592 10(46) x 10 -10 (±3.9 ppm), which agrees with 
the standard model to within one standard deviation. These results are displayed 
graphically in Fig. || Also shown is the projected error from the 1999 data, and the 
±0.35 ppm goal of E821. 



5 Outlook and Conclusions 



The standard model has been remarkably successful in describing a wide range of 
phenomena including the entire body of results from LEP. The new result from E821, 
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Measurement 


Value xlO 10 


CERN|IJ (j+ 
CERN[|ig] ir 
E821 133| //+ 7r inj 
E821* fi+ /i inj 


116 591 03 (120) (10 ppm) 
116 593 65 (120) (10 ppm) 
116 592 51 (150) (13 ppm) 
116 591 91 (59) ( 5 ppm) 


New World Average 


116 592 10 (46) ( 4 ppm) 



Table 1: The values of from the four most precise experiments. The new value 
for n^l fi p has been used to get a M from the measured ratio uj a /uj p . For the average, 
X 2 /v — 0.92. The goal of the experiment is an error of ±4.0 in the units above. 
* Preliminary 



116 595 000- - 



^ 116 594 000- - 



x 116 593 000--° 



116 592 000- 
Theory 

116 591 000- 



116 590 000-- 



E 

Oh 
Oh 



Oh 



— <) 



Oh 
Oh 



Oh 
Oh 



&H 



+ 1 

! 





U EC 



■ w 

■ u 



00 

w 



C-4 
CO 

pq 



0) 

bO 
ffl 

HH 

0) 

> 

0) 



n3 

HHJ 

-O 

ON 
ON 

CTn 



* preliminary 
*projected errors with current value 



Figure 8: The four precise measurements of and the standard model value. To 
show the potential improvement available, the projected errors of 1 ppm and 0.35 
ppm are shown on the right hand side with the central value equal to the present 
world average. 



which has a precision of 5 parts per million, lowers the uncertainty on our knowledge 
of (g — 2)n to 3.9 ppm. At that level there is good agreement with the standard model, 
but there is still room for the observation of a value at the level of one ppm or better 
which could agree with the previous measurements, and disagree with the standard 
model. We expect our one ppm result to be available before the end of 2000. 



9 



This experiment is very much a work in progress. The initial design of E821 was 
made with the systematic error budget of 0.12 ppm. Many of the improvements which 
were made to the CERN technique have indeed worked, making it straightforward to 
obtain a systematic error less than one ppm. While we have been learning about our 
systematic errors from the beginning of the first pion injection run, with the sub-ppm 
statistics available with muon injection we are now challenged to push the analysis 
of systematic errors to the limit. Our ultimate goal is a statistical answer at the 0.3 
ppm level, with systematic errors at perhaps half this level. Our work thus far seems 
to indicate that this level of systematic error will be possible. 

A new inflector magnet has been installed in the ring, which has a fringe field in 
the storage region a factor of 5 less than the old inflector had. This will improve our 
ability to map the field everywhere and will help to reduce the uncertainty on our 
knowledge of the field. The principal magnetic field issue remaining is our ability to 
track the field with time, and to understand the full calibration of the NMR probes 
to a few tenths of a ppm. 

The other principal challenge is how to handle pile-up in the detectors without 
discarding an unreasonable portion of the data set. The dominant analysis effort now 
is to understand this effect. 

It has been 20 years since the CERN experiment presented its final report which 
verified the hadronic contribution at the eight standard deviation level, but which left 
it to a future experiment to verify the electroweak contribution. The new Brookhaven 
experiment is now approaching that goal. Within the next few years either the 
standard model value will be confirmed, or evidence of a new contribution to the 
muon (g — 2) will be discovered. 

I wish to acknowledge the efforts of the many collaborators who have worked on 
(g — 2) over the past 15 years. The steady support of the funding agencies (US and 
abroad) and the Brookhaven Laboratory management was essential to our reaching 
this point. I wish to thank R. Carey, D. Hertzog, K. Jungmann, V. Hughes, J. Miller, 
Y. Semertzidis and E. Sichtermann, for their comments on this manuscript. 
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